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ABSTRACT 

We study the Hi A-band Tully-Fisher relation and the baryonic Tully-Fisher relation for a sample of 16 early-type galaxies, taken from 
the ATLAS^*^ sample, which all have very regular Hi disks extending well beyond the optical body (> 5 Rcs)- We use the kinematics 
of these disks to estimate the circular velocity at large radii for these galaxies. We find that the Tully-Fisher relation for our early-type 
galaxies is offset by about 0.5-0.7 mag from the relation for spiral galaxies, in the sense that early-type galaxies are dimmer for a given 
circular velocity. The residuals with respect to the spiral Tully-Fisher relation correlate with estimates of the stellar mass-to-light ratio, 
suggesting that the offset between the relations is mainly driven by differences in stellar populations. 

We also observe a small offset between our Tully-Fisher relation with the relation derived for the ATLASsample based on CO data 
representing the galaxies’ inner regions (< 1 Reir)- This indicates that the circular velocities at large radii are systematically 10% lower 
than those near 0.5-1 Rcs, in fine with recent determinations of the shape of the mass profile of early-type galaxies. 

The baryonic Tully-Fisher relation of our sample is distinctly tighter than the standard one, in particular when using mass-to-light 
ratios based on dynamical models of the stellar kinematics. We find that the early-type galaxies fall on the spiral baryonic Tully-Fisher 
relation if one assumes MJLk = 0.54 MqILq for the stellar populations of the spirals, a value similar to that found by recent studies 
of the dynamics of spiral galaxies. Such a mass-to-light ratio for spiral galaxies would imply that their disks are 60-70% of maximal. 
Our analysis increases the range of galaxy morphologies for which the baryonic Tully-Fisher relations holds, strengthening previous 
claims that it is a more fundamental scaling relation than the classical Tully-Fisher relation. 

Key words, galaxies: kinematics and dynamics - galaxies: elliptical and lenticular 


1. Introduction 

The Tully-Fisher relation (hereafter, TER; iTullv & FisherllTo?^ 
is an established fundamental scaling relation for spiral galax¬ 
ies. It was originally discovered as a tight correlation between 
the width of the integrated H i spectrum and the absolute mag¬ 
nitude, but is primarily an empirical relation between circular 
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velocity and luminosity. Explaining the TER is fundamental for 
understanding how galaxies form and evolve. Moreover, since 
luminosity is distance dependent, while circular velocity is not, 
a historically important application of the TER is the determina¬ 
tion of distances in the local Universe. 

Since its discovery, there has been significant debate about 
how to interpret the TER in terms of more fundamental galaxy 
properties such as total stellar mass. For this reason, modern 
TER studies tend to observe galaxies at infrared wavelengths, 
where luminosity is a better tracer of stellar mass. An even 
more fundamental relation might exist between the circular ve- 
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locity and total baryon ic mass, the baryonic TFR (IWalkejll999t 
iMcGaugh et al.l I 2 OOOI) . For example, while the TFR becomes 
less-well de fined for smaller g alaxies (circular yelocities below 
50 km s '), iMcGaugtil (l2012li shows that these systems follow 
the same baryonic TFR as brighter ones. 

The study of the TFR as a function of galaxy morphological 
type can clarify to what extent the slope, intercept and scatter 
of the relation are determined by stellar MjL or shape of the 
mass distribution. Earlier studies indicate that early-type spi¬ 
rals are fainter for a giyen circular yelocity than la te-type spi¬ 
rals (lRobertslll978t iRusselll 120041 IShen et al.ll200'^ . Howeyer, 
this difference reduces when stellar population effects are min¬ 
imised by considering near-infrared instead of optical photom¬ 
etry ( e.g., lAaronson & MouldlllOSSt Noordermeer & VerheiienI 
l2007h . One other important finding is that the shape of the TFR, 
in particular whether the TFR exhibits a ’kink’ towards high 
circular yelocities for more luminous galaxies, depends on how 
and at which radius the circular yelocities are measured, reflect¬ 
ing systematic differences in t he shape of the mass distribution 
and hence of the rotation curye ( Noo rdermeer & Verheiienl200'^ 


Williams, Bureau & Cappellari 201(iir 


Recently, TFR studies haye been extended to early-type 
galaxies (E and SO, hereafter ETGs). These objects are impor¬ 
tant because by comparing their TER to that of spirals we haye 
a larger ’baseline’ for inyestigating how the stellar mass-to-light 
ratio and the shape of the total mass distribution influ ence the 
TER. iBedregal, Aragon-Salamanca & Merrifieldl (l2006ll studied 
the B- and /T-band TER of lenticular (SO) galaxies and found 
that in both bands the TER for SO’s lies below the TER for spiral 
galaxies (1.2 mag i n jiT-band). Hence, SO galaxies r otate faster 
and/or are dimmer. IWilliams, Bureau & Cappellari! (1201 Oh also 
found an offset based on a study of 14 spirals and 14 SOs using 
yario us yelocity tracers an d estimate it to be 0.53 mag in the Ks- 
band. l^rtesi et al] (12013h performed a study of the SO TER us¬ 
ing stellar kinematics from planetary nebulae spectra and found 
the SO TER to be offs et from the spiral T ER by a full magnitude 
in ^f-band. Similarly. [Rawle et al.l (1201 3h find an offset in the 
band between the TER of sp irals and ETGs in the Coma cluster 
of about 0.8 mag. Recently. iJaffe et ^ (12014h used the ionised 
gas emission of 24 ETGs to construct a TER and found that their 
ETGs are 1.7 mag fainter than spirals in the B-band. 

Until recently, studies of the ETG TER haye been challeng¬ 
ing because such galaxies are typically gas-poor compared to 
spirals and measuring their circular yelocity is less straightfor¬ 
ward, in particular at large radii. Howeyer, recent suryeys, no¬ 
tably the SAURON (Ide Zeeuw et alJ l2002h and the ATLAS® 
projects dCappellari et al.l 2011 ). haye shown that a significant 


fraction of ETGs host significant amounts of molecular and 


or ring (iMorganti et all 

2006HOoster 

00 et alJl2010t Young et alJ 

120111 ISerra etalJ 2012 

; Dayis et al.l 

20131; lAlatalo et al.ll2013li. 

As part of the ATLAS^*^ project. Dayis et al. (|201 l|i studied the 


iMorganti et al.l (l2006h . lOosterloo et alJ (l2010h and ISerra et all 
(1201 2h around ETGs, and study the TER of these galaxies us¬ 
ing circular yelocities measured at yery large radius (> 5Reff)- 
This paper is structured as follows. In Section 2 we describe 
the sample and Hi data. In Section 3 we discuss the methods 
used to estimate the H i circular yelocity. In Section 4 we present 
the HI jT-band TE R and a comparison with the CO TER from 
iDayis et aP (1201 Ih and with TERs for spiral galaxies. In Section 
5 we present the Hi baryonic TER and in Section 6 we sum¬ 
marise our findings and present the conclusions. 

2. Sample and data 

This work is based on the data from the ATLAS^*^ Hi sur- 
yey of ETGs carried out with the Westerbo rk Synthesis Ra¬ 
dio Telescope (WSRT) by ISerra et ^ (l2012h . This suryey tar¬ 
geted a yolume-limited sample of 166 morphologically selected 
ETGs with distance D < 42 Mpc and brighter than Mk - 
-21.5. This sample consists of those galaxies of the larger, 
yolume-limited ATLAS^*^ sample w hich are obseryable wit h the 
WSRT and includes the Hi data of Morgantie^alJ (2006) and 
Oosterloo et al.l (2010h . We refer to ICappellari et alJ ( 2011 ) and 


Serra et al.l (12012h for more details on the sample selection and 


CO TER and their most important findings are that in many 
ETGs CO can be used as a tracer for the circular yelocity of 
the flat part of the rotation curye and that the CO TER is off¬ 
set from the TER for spirals by abo ut 1 mag in /T-band. The 
CO disks used by iDayis et ^ (1201 ll) typically extend to about 
1 Reff and, therefore, the TER in their work reflects the situ- 
atio n for the inner regions of ETGs . As has been found by, 
e.g.. I Noordermeer & VerheiienI (1200'^ . the shape of the obseryed 
TER depends on at which radius the circular yelocities are mea¬ 
sured, reflecting the fact that the shape of rotation curyes yaries 
systematically with galaxy mass and morphology. Here we take 
adyantage of the large size of many of the H i discs found by 


on the Hi obse r yation s. 

ISerra et ^ (1201 2h report the detection of H i in 53 ETGs. In 
particular, in 34 objects the H i is found to form a rotating disc or 
ring with radii yarying from a few to tens of kpc (classes d and D 
in that paper). These systems exh ibit relatiyely regular H i kine¬ 
matics (see also lSerra et al.l20T4l for a discussion of the kinemat¬ 
ics of these H i disks) and are therefore the starting point for our 
TER study as these H i disks offer the possibility to deriye accu¬ 
rate rotation yelocities. A second step in the data selection has 
been to perform a harmonic decomposition of the yelocity field 
of the oute r regions of these 34 g alaxies using the KINEMETRY 
software of lKraino yic et al. l (l2006h. T he yelocity fields used are 
those published by ISerra et alJ (l2014l) . A harmonic decomposi¬ 
tion of the yelocity field giyes information about to what extent 
circular rotation dominates the obseryed kinematics compared 
to other, non-circular motions. Our selection criterion was that 
the amplitude of the non-circular motions must be less than 10% 
of the circular rotation. This eliminates cases where the H i in 
the outer regions is not sufficiently settled in a disk, or galaxies 
where strong streaming motions due to a bar might affect the re¬ 
sults. Moreoyer, since we aim to study the TER using yelocities 
obtained at large radius, we exclude galaxies in which the Hi 
component is not yery extended compared to the optical size. 

These criteria reduce the sample of 34 objects by exclud¬ 
ing galaxies with signs of interaction (NGC 3619, NGC 5103, 
NGC 5173), galaxies with complex Hi kinematics (e.g., 
strong warps or non-circular motions; NGC 2594, NGC 2764, 
NGC 4036, NGC 5631, UGC 03960, UGC 09519) and galax¬ 
ies whose Hi disc is smaller than or comparable to the stel¬ 
lar body (NGC 3032, NGC 3182, NGC 3414, NGC 3489, 
NGC 3499, NGC 4150, NGC 4710, NGC 5422, NGC 5866, 
UGC 05408). This leayes us with 15 galaxies. We add to this 
sample NGC 2974, whi ch is part of the A TLAS^*^ sample but 
not of the Hi suryey of ISerra et al.l (1201 2h because of its low 
declination. 

Since we consider only relatiyely H i-rich galaxies with reg¬ 
ular kinematics, we inyestigate whether this introduces any bias. 
In Eig. [H we show the mass-size relation of the galaxies from 
ISerra et ^ (1201 2h with the sub-sample selected by us for thi s 
paper indicated (the data are taken from ICappellari et~^l2013ll . 
One can see that yery massiye and yery small galaxies appear 
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logic AfjAM (Mq) 

Fig. 1. Mass-size plane for the galaxies from the s ample in this pape r 
(the red symbols) compared to the full sample from ISerra et alj (|2012|) . 
Data are taken from Cappellari et al. (2013b) while the definitio n of 
the HI classes is very similar to the one used bv ISerra et al.l ( l2012h : D: 
large H i disk (larger than optical body), R large H i ring (larger than op¬ 
tical body), d'. small H i disk (confined within optical body), n: unsettled 
(tails, streams), c: scattered Hi clouds. The blue, dashed lines give loci 
of constant velocity dispersion. 


to be missing from our sub-sample. The main consequence is 
that the range of circular velocities over which we can study the 
TFR is somewhat limited. Furthermore, Cappellari et al. (2013a) 
showed that the bulge-to-disc ratio, which affects the shape of 
the rotation curve, changes systematically with velocity disper¬ 
sion (represented by the blue, dashed lines in Fig. 1). The figure 
shows that our TFR sample nearly covers the entire range of ve¬ 
locity dispersion of the full sample. Therefore, we conclude that 
our selection has not introduced any major bias. 

In order to construct a TFR, we need the H i circular veloc¬ 
ities for these 1 6 galaxies. F or NGC 2685 and NGC 2974 we 
use the results of l.T6zsa et al.l (120091) and IWeii mans et aP (l2008l) . 
respectively. For the remaining 14 galaxies we use the methods 
described in Section 3. Since the finite spatial resolution of our 
observations gives rise to beam-smearing in the inner regions of 
the velocity fields, we do not derive complete rotation curves for 
the galaxies. Instead, we focus on measuring the circular veloc¬ 
ity at large radius since this is the relevant quantity for our TFR. 

For the sample ga laxies we adopt t he abs olute magnitudes 
and distances listed in ICappellari et al.] (1201 lb . We refer to that 
paper for details on the photometry and on the distance estimate 
of all galaxies. We assume that the main uncertainties in the K- 
band magnitudes, and of the stellar masses we derive from these 
magnitudes, are caused by distance uncertainties. We estimate 
these from the average deviation of redshift distances from the 
NED-D redshift-independent distance^ which we find to be ap¬ 
proximately 21 percent. This implies an uncertainty in the abso¬ 
lute magnitudes of 0.41 mag. 


3. Methods 

An important aspect of the TFR is what quantity should be 
used as velocity. Historically, the TFR has been constructed 
using the width of the integrated Hi spectrum. However, 
more recent studies tend to use directly the circular velocity, 

* http://nedwww.ipac.caltech.edu/Library/Distances/ 


which can be estimated using a variety of metho ds includ¬ 
ing H I, CO and Hg resolved rotation curves (e.g., [Verheii^ 
|2^ iDavis et all 120111 i McGaugh Rubin & de Blo3 ~ l20m 
planetary nebulae (|Cortesi_eLalJ l2013l) or dynamical modelling 
iwilliams. Bureau & Cappellarii200^ . These different methods 
may result in estimates of the circular yelocity at different points 
along the rotation curye, and it is important to check that inde¬ 
pendent methods deliyer consistent results. Recent analyses haye 
indeed highli ghted the importance of the rotati on curye shape for 
TFR studies (INoordermeer & Verheiienll2007h . We reiterate that 
in this paper we focus on circular yelocities for radii much larger 
than the optical size of the galaxies. 

To estimate the H i circular yelocity, we start by analysing the 
moment-1 yelocity fie lds produced by th e WSRT data reduction 
pipeline described in ISerra et ^ (1201 2h and that are shown in 
ISerra et ^ (l2014l) . Howeyer, these yelocity fields do not always 
represent the H i kinematics accurately, in particular for objects 
with a projected size small enough to suffer from beam-smearing 
or for galaxies where the outer disk is significantly warped. For 
these objects we attempt to take beam-smearing effects into ac¬ 
count by analysing the H i data cube directly. Which technique 
we haye used for which galaxy is indicated in Table [1] We dis¬ 
cuss the analysis procedure applied to yelocity fields and to H i 
cubes in Sections [3.11 and [T2l respectiyely. Table [T] lists the re¬ 
sulting circular yelocity of all galaxies in our sample. As men¬ 
tioned in Section|2] giyen the size of the H i beam relatiye to the 
extent of the galaxies, we do not attempt to deriye complete rota¬ 
tion curyes coyering all radii. To construct the TFR, we focus on 
denying a single circular yelocity for the outermost point of the 
rotation curye. The radius at which we were able to determine 
this yelocity differs from galaxy to galaxy and ranges from 8 to 
28 kpc (15 kpc on ayerage) or R/Reg from 3.4 to 13.7 (7.3 on 
ayerage). 


3.1. Velocity field analysis 

For galaxies for which the H i is well resolyed, we deriye the H i 
circular yelocity at large radius b y fitting tilted-ring models to 
the yelocity fiel d (lBegemanlll987h using t he rotcur task in the 
GIPSY package (|yan der Hulst et al]ll992h . This task performs a 
least-square fit to the yelocity field by modelling the galaxy as 
a set of rings with increasing radius while yarying the following 
parameters of each ring: centre, systemic yelocity, inclination, 
position angle and circular yelocity. We adopt a spacing between 
the rings equal to the H i beam size. As starting parameters for 
the fitting, the inclination is based on the axis ratio of the Hi 
column density map, while the initial position angle, circular- 
and systemic yelocities are estimated by using position-yelocity 
diagrams. In the fitting we used the following strategy. Firstly, 
we yary only the coordinates of the centre for all rings together 
and then fix them for the rest of the fitting procedure. Secondly, 
we fit all parameters and fix the systemic yelocity of subsequent 
fits to one yalue for all rings based on this fit. Finally, we solye 
for the inclination, position angle and circular yelocity (Wcirc) for 
the outer rings. 

The uncertainty in the circular yelocity is largely determined 
by the error in the inclination. It is therefore important to esti¬ 
mate the latter. To be able to assess the accuracy of the best¬ 
fitting parameters we follow a two-step strategy. First, we fit, for 
each ring, higher order harmonic terms to the yelocity field: 

v(R,9) - tisys + "^CkiR) cos kO + Sk(R)smk0, (1) 

k 
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Fig. 2. Left: NGC 3941 column density map. The contours are 0.5, 1, 2, 4 x 10^° cm“^. The ellipse on the bottom-left describes the the beam 
shape. Right: First-moment velocity field. The contours range from 780 to 1080 km/s with 30 km/s steps. The overlaid ellipse describes the model 
geometry (kinematic centre, position angle, inclination and radius). The ellipse on the bottom-left describes the the beam shape. 


where v^sys is the systemic velocity, Ck(R) and SkiR) are 
the harmonic amplitudes of order k and 6 is the azimuthal 
angle in the plane of the galaxy and Ccirc = ci / sin i. 
We use the first four orders of this expansion. This fit 
is perfor med using the routine reswri from the GIPSY 
package d Schoenmakers, Fra nx & de Zeeuwl Il997h . Following 
iFranx. vm Gorkom & de Zeeuwfil 19941) . we use the fact that an 
error in the inclination angle can be directly related to the C 3 - 
term. We make use of this by estimating the error in the inclina¬ 
tion by constructing velocity fields for the outer rings for a range 
of inclinations around the value obtained from the tilted-ring fit 
and see for which inclinations we detect a non-zero C 3 term in 
the difference between the observed and model velocity field. 

Summarising our method, we first fit a regular tilted-ring 
model to estimate the rings’ centre, systemic velocity, position 
angle and inclination as described above. We then use these val¬ 
ues as input for the harmonic decomposition, which we run mul¬ 
tiple times on the residual velocity fields that we derive by vary¬ 
ing each time the inclination of the corresponding model velocity 
field in steps of 1 deg. From these harmonic decompositions, we 
derive the function cj,(i) to determine the uncertainty on inclina¬ 
tion. 

We illustrate this technique showing its application to the ve¬ 
locity field of NGC 3941. This galaxy hosts an H i ring which is 
well resolved, as shown in Fig.|2 The best-fitting parameters re¬ 
sulting from a tilted-ring fit of the outer rings are i - 57 deg and 
i^drc = 148 km s“*. Figure [3 shows the value of C 3 as a func¬ 
tion of inclination offset from the best fitting inclination of 57 
deg and has, of course, its minimum at the best-fitting inclina¬ 
tion. The task reswri returns an error bar 0-3 on C3. We define 
the uncertainty on i as the interval Ai within which C 3 is equal to 
its minimum value within one 0 - 3 . The horizontal line in Fig. [3 
represents this definition and, in this case, returns an error on i 
of ± 1 deg. 

We can now use this result to estimate the uncertainty on 
the circular velocity. In this case, running rotcur again at fixed 
inclinations 56 deg and 58 deg respectively, and solving for the 



A* (°) 


Fig. 3. NGC 3941 inclination angle offset versus C 3 -amplitude. The hor¬ 
izontal line describes the csfbestfit) - 1 - Acsfbestfit), where Acs is the l-cr 
uncertainty of the Cs-term. Since the minimum occurs at the best-fitting 
inclination angle, and any offset from the best-fitting inclination gives 
rise to a C 3 amplitude, the inclination can be confidently constrained. 


circular velocity, yields a change of +3 km s ' on the latter. In 
fact, in cases like this we prefer to assign a conservative estimate 
of +8 km s * to the Ccirc uncertainty, corresponding to about half 
the velocity resolution. 

3.2. Model data cubes 

For poorly-resolved galaxies, or galaxies with a significantly 
warped disk, beam-smearing persists even in the outer regions 
and we model the full H i data cube instead of the velocity field 
(NGC 2824, NGC 3838, NGC 4203). In these cases we there¬ 
fore perform a tilted-ring analysis using the TiRiFiC software 
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Fig. 4. Left. NGC 3838 column density map. The contours are 5, 10 x 10** cm“^. The ellipse on the bottom-right describes the the beam shape. 
Right. First-moment velocity field. The contours range from 1100 to 1450 km/s with 20 km/s steps. The ellipse on the bottom-left describes the 
beam shape. 


jJozsa et alJl2007lPI which is specifically designed to take beam 
smearing effects into account. Our goal is not to construct a 
perfectly matching kinematical model of the entire galaxy, but 
rather to construct simple models to constrain the inclination 
and circular velocity of the galaxy for the outermost points of 
the rotation curve. Since the TiRiFiC opitimiser is sensitive to 
local minima, accurate initial guess values are required for 
the parameters. The starting values for the fit are chosen analo¬ 
gous to Section 13.11 After setting the initial parameters, we fix 
some parameters and fitting only a few parameters at a time. This 
procedure is done iteratively until a stable solution emerges. For 
simplicity, and considering the limited scope of our fitting (i.e., 
focus on the outermost point of the rotation curve), we fit po¬ 
sition angle and inclination as being constant with radius if no 
strong warp is present. In the case of NGC 4203, however, a 
warp is present and this warp is modelled as a linear function in 
inclination and position angle. 

Although the velocity fields do not allow for a direct deriva¬ 
tion of the tilted-ring parameters for poorly resolved galaxies, 
we can still use them to estimate the uncertainties on the incli¬ 
nation angle derived from the data cube and therefore the un¬ 
certainty on the circular velocity. In practice, the uncertainty on 
^circ is constrained by deriving a set of model velocity fields us¬ 
ing the GIPSY-task velfi with inclinations offset from the best¬ 
fitting inclination angle. By inspecting the iso-velocity contours 
of these velocity fields, we estimate a confidence interval for the 
inclination and from this we estimate the uncertainty on Vdrc- 

We illustrate this method showing its application to NGC 
3838. The Hi in NGC 3838 is faint and not well resolved by 
the WSRT beam, causing a large uncertainty in the intensity- 
weighted velocity. Fig, shows the column density map. The 
ring as reported bv ISerra et al.l (1201 2h is not resolved along the 
minor axis and therefore a basic rotcur tilted-ring fit and a har¬ 
monic decomposition of the velocity field will not yield any re¬ 
liable results. Using fits to the data cube, we find a best fitting 


^ http://gigjozsa.github.io/tirific/ 



Angular Offset 


Fig. 5. NGC 3838 position-velocity diagram. The blue contour corre¬ 
sponds to the observation and the red contour to the model. In both 
cases, the countour level is at 0.8 mjy beam^*, which is approximately 
the 2 - 0 - noise level of the data cube. 


inclination and circular velocity of 66 deg and 159 km s“’, re¬ 
spectively. We derive the uncertainty on the latter value by build¬ 
ing a number of H i velocity fields for a range of inclination val¬ 
ues around the best-fitting one (changing inclination in steps of 
1 deg). We estimate by eye the acceptable range of inclination, 
which in this case is +8 deg. Finally, we again generate a model 
data cube having fixed the inclination at the boundaries of this 
interval and obtain an uncertainty on the circular velocity of +16 
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Table 1. The data relevant for the TFRs discussed in this paper. 


ID 

(1) 

D 

(Mpc) 

(2) 

Method 

(3) 

i 

(deg.) 

(4) 

Ai 

(deg.) 

(5) 

^circ 

(kms^') 

(6) 

AUcirc 

(kms"^) 

(7) 

log Mhi 
(M o) 

(8) 

logLr 

L© 

(9) 

log 

MqI 

(10) 

logM/4^™ 

Mo/Lor 

(11) 

NGC 2685 

16.7 

c 

65 

4 

144 

10 

9.33 

9.857 

0.455 

0.620 

NGC 2824 

40.7 

c 

65 

8 

162 

10 

7.59 

9.889 

0.628 

0.343 

NGC 2859 

27 

V 

35 

5 

215 

41 

8.46 

10.404 

0.568 

0.732 

NGC 2974 

20.9 

V 

60 

2 

310 

10 

8.74 

10.152 

0.981 

0.787 

NGC 3522 

25.5 

V 

41 

2 

121 

8 

8.47 

9.600 

0.705 

0.648 

NGC 3626 

19.5 

V 

54 

1 

169 

8 

8.94 

10.102 

0.437 

0.261 

NGC 3838 

23.5 

c 

66 

8 

159 

14 

8.38 

9.772 

0.589 

0.694 

NGC 3941 

11.9 

V 

57 

1 

148 

8 

8.73 

9.940 

0.400 

0.693 

NGC 3945 

23.2 

V 

57 

-5,-bl 

237 

13 

8.85 

10.394 

0.628 

0.760 

NGC 3998 

13.7 

V 

66 

2 

246 

20 

8.45 

9.967 

0.971 

0.818 

NGC 4203 

14.7 

c 

30 

3 

197 

35 

9.15 

10.067 

0.537 

0.827 

NGC 4262 

15.4 

V 

67 

-8 ,-hI 

198 

10 

8.69 

9.726 

0.753 

0.782 

NGC 4278 

15.6 

V 

45 

-2,-h5 

256 

26 

8.80 

10.247 

0.829 

0.846 

NGC 5582 

27.7 

V 

50 

2 

258 

10 

9.65 

10.140 

0.722 

0.752 

NGC 6798 

37.5 

V 

52 

2 

190 

8 

9.38 

10.028 

0.660 

0.684 

UGC6176 

40.1 

V 

47 

-2,+6 

144 

14 

9.02 

9.751 

0.685 

0.328 


Notes. List of the galaxies where the circular velocit ies at large radii are de rived from the Hi kinematics. Column (1): the name is the galaxy 
designation. Column (2); galaxy distance as listed in ICappellari et all i201lh . Column (3): c indicates that the circular velocity is based on an 
analysis of the full data cube, whereas v indicates a velocity-field based circular velo city. Column (4) a nd (5): inclination angle and uncertainty. 
Column (6) and (7): circular velocity and uncertainty. Column (8): Hi-mass ffrom ISerra et a J j2012|)). Column (9) : r-band luminosity (from 
ICappellari et alj|2013ll . Column (10): r-band mass-to-light ratios from stellar dynamics (from Cappellari et alj|2013l) and Column (11): r-band 
mass-to-light ratios from star formation histories (from Cappellari et al. 2013a). 



Fig. 6. The H i Tully-Fisher relation for our 16 galaxies (red circles; W = 
Zurire)- The black sy mbols give the data based on pla netary nebulae from 
ICortesi et alJ i2013h for comparison. The CO TFR of|Da^s_et^^j20nh 
is give n by the drawn line. The dashed line gives the iTullv&Pierea 
ll2000h /f-band TFR for spiral galaxies. 


km s Figure |5] shows a comparison of the observed position- 
velocity diagram with the one based on the model. 


4. The HI K-band Tully-Fisher relation of early-type 
galaxies 

In this section we discuss the /T-band TFR o f ETGs, compare it 
to the one derived using CO kinematics by iDavis et ^ (1201 ih 
and compare our results to the TFR of spiral galaxies. The rele¬ 
vant data used are summarised in Table [1] 

Figure|6]shows the H i TFR as found from our data. For com¬ 
parison, the hgure includes the early-type galaxies whose circu¬ 
lar vel ocity is estimated using planetary nebulae bv ICortesi et al.l 
(l2013l) . Although we do not use these galaxies in the formal TFR 
fits below, the hgure shows that the two samples are fully con¬ 
sistent. 

We determine the parameters of the TFR by htting a relation 
of the form 

Mk = a{logW-2.6) + b (2) 

to our data. We perform such hts a number of times 
with different constraints using a weighted least-square ht 
of the inverse relation usin g the MPFITEXY package of 
I Williams. Bureau & Cappella rii (1201 Oh . which depends on the 
MPFIT-nackage (lMarkwardlil2009l) . We list the results in Table 
12] Table 121 shows that the scatter around the unconstrained ht is 
only slightly larger than expected from errors in the data (mainly 
distance errors). We also hnd that the slope a and intercept b 
are not well constrained by such a ht. This is because the cor¬ 
relation p between the points is only modest - the un-weighted 
correlation coefficient is 0.74. Therefore, to be able to make a 
better comparison with other published TFRs, we have also per¬ 
formed hts where we have constrained the slope to that of other 
published TFRs. 
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Table 2. Parameters of the H i K-band and baryonic TFR. 


TFR 

slope a 

intercept b 

cr 

P 

A'-band 

unconstrained 

-4.67 + 0.99 

-23.25 + 0.12 

0.43 

0.74 

slope fixed 

-8.38 (CO) 

-23.31 + 0.17 

0.67 


slope fixed 

-8.78 (TPOO) 

-23.32 + 0.18 

0.69 


slope fixed 

-9.13 (WBCIO) 

-23.32 + 0.19 

0.72 


slope fixed 

-9.64 (NV07) 

-23.33 + 0.20 

0.77 


Baryonic 

unconstrained (JAM) 

2.38+0.25 

10.71+0.03 

0.24 

0.94 

unconstrained (SFH) 

2.51 +0.42 

10.71+0.05 

0.47 

0.80 

slope fixed (JAM) 

3.36 (NV07) 

10.72 + 0.04 

0.39 


slope fixed (SFH) 

3.36 (NV07) 

10.72 + 0.06 

0.57 


slope fixed (JAM) 

3.84 (McG12) 

10.73 + 0.05 

0.50 



Notes. The rows in the table present the best-fitting parameters for the A'-band TFR (all inverse fits) and the JAM baryonic TF R for the 16 galaxie s 
in Tab le [H The slopes withou t uncertain ties correspond to the cases where the slope ha s been fixed to a reference TFR: C O: IPavis et al. j[20nL 
TPOO: iTullv & Pierce! dTOOOh . WBCIO: I Williams. Bureau & Cappellaril ( 1201 Oh . NV07: iNoordermeer & VerheiienI ( l2007h and McG: iMcGaughl 
(l2012h . The scatter cr in mass of the baryonic TFR has been converted into magnitudes by multiplying by 2.5. The last column gives the correlation 
coefficient of the data. 


4.1. Comparison to the CO K-band TFR 

We start by comparing t he H i K-band TFR of ETGs to the anal¬ 
ogous CO relation from iDavis et al.l (l201 ih . This may be inter¬ 
esting because, since the typical CO disk extends only to about 
1 /?eff while our H i data extends much further, it allows to com¬ 
pare the TFR valid for the inner regions of ETGs with the one of 
the entire galaxy. This, in turn, m ay tell us somethin g about the 
shape of the mass distribution. In iDavis et aP (1201 Ih . a number 
of TFR relations are derived using different samples and meth¬ 
ods. Here we compare our data to their “hybrid’ TFR, which is 
based on the largest sample of galaxies for which either single¬ 
dish or interferometric CO data suggest that the molecular gas 
has reached the flat part of the rotation curve. 

We investigate a potential offset between H i and CO TFR by 
keeping the slope a fixed to the CO value of -8.38. The resulting 
intercept isb — -23.31 +0.17 (Table|2]i, to be compared with the 
CO value of -23.12 ± 0.09. Therefore, the difference between 
CO- and H i TFR intercept is constrained to be between +0.45 
and -0.07. At fixed M^, this corresponds to a ratio of H i to CO 
circular velocity between 0.88 and 1.02. Therefore, the circular 
velocity of a typical ETG decreases between 0 and 10 percent 
when going from 0.5 - 1 Reft, where it is traced by CO , to the 
outer regions, where it is traced by H i. ICapnellari et al.l (l2015h 
recently found that the mass profiles of a sample of 14 ETGs can 
be well described by a nearly isothermal power law ptot k 
with an average logarithmic slope y = 2.19 with remarkable 
small scatter in y. The small decrease in circular velocity we 
observe is consistent with such mass profiles. We will explore the 
mass profiles of our sample galaxies, combining the ATFAS^° 
data on the stellar kinematics with our H i data, in more detail in 
a future paper (Serra et al. in prep.). 

4.2. Comparison to the K-band TFR for spirais 

Another interesting exercise i s to compare our d ata with the TFR 
found for later-type galaxies. IPavis et ^ (1201 ih reported an off¬ 
set between their CO TFR for ETGs of the ATFAS^° sample 
with the one for spiral galaxies in the sense that for a given cir¬ 
cular velocity, ETGs are fainter. Our data also indicate such an 
offset (Fig. @1. However, since our H i TFR for early-type galax¬ 
ies is slightly offset from the one based on CO data in the other 



log(W) (km s ^) 


Fig- 7- Our data (red symbols) compared with the TFR derived by 
INoordermeer & VerheiienI ( 120071) using the asymptotic circular veloc¬ 
ities (black symbols). 


direction, the offset from the spiral relation is somewhat smaller. 
If we fit a relation wit h the slope fixed to tha t of the TFR found 
for spiral galaxies by iTullv & Piercel (l2000l) . we find an offset 
from their relat ion of 0.70 magni tudes, compared to 0.93 mag 
found found bv iDavis et al.l (1201 ll) . We therefore find a shift be¬ 
tween the relations of about 0.2 mag. 

We can also compare our r esults with those of 
IWilliams. Bureau & Cappellaril (1201 Oh . who studied an en¬ 
semble of 14 SO and 14 spiral galaxies and found their SO TFR 
to be 0.5 mag below thei r spiral TFR. A fit to our data fix i ng the 
slope to that found by IWilliams. Bureau & Cappellaril (1201 Oh 
gives a 0.25 mag offset from their sample of spirals. 
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Fi g. 8 . Residual (in magnitude) of ou r galaxies from the TFR published 
bv lNoordermeer & Verheij^ (|2007| ) as fu nction of mass-to-light ratio 
(r band) derived bv ICaPDellari et all ll2013li using dynamical modelling. 


Finally, iNoordermeer & VerheiienI (l2007h study the /f-band 
TFR of massive spirals. Noordermeer & Verheijen find a ^f-band 
TFR of the form Mk = (-9.64 + 0.25)(log IT - 2.6) + (-23.95 + 
0.04). A fit to our data with the sl ope fixed to the slope found by 
Noordermeer & VerheiienI (l2007h (see Table |2]i shows an offset 
of 0.62 mag. 

Our analysis appears to confirm a zero-point offset of 
the /T-band ETG TFR with respec t to th e TFR for spirals 
of about 0.5-0.7 mag. iDavis et 51 ( 2011h sugges ted that a 


simple fading mechanism (see e.g. iDresslej Il980l) combined 
with low levels of residual star formation (as implied by 
the presence of molecular gas ) woul d explain such an offset. 
Williams. Bureau & Cappellaril (l2010li remark that a difference 
in size between ETGs and spirals could provide an alternative 
explanations for the offset. Their work showed that an offset be¬ 
tween their SO TER and spiral TER is also present in the dynam¬ 
ical mass TER. Dynamical mass scales as where r is some 
characteristic radius. Therefore, if ETGs were more concentrated 
than spirals, an offset would occur between the relations for these 
two types of galaxies. Below we show that no significant offset 
is present in the baryonic TER, suggesting that the differences in 
mass distribution do not play an important role. 

In interpreting these numbers, one should take into ac¬ 
count that the various studies use velocities measured in dif¬ 
ferent wa ys and for different regions of the ga laxies. Eor ex¬ 
ample, the lWilliams, Bureau & Cappell aril ( 20 lOl) analysis refers 
to the cental regions while that of iNoordermeer & VerheiienI 
(l2()0^ uses F[ I circular velocities at large radius and is more 
similar to ours. In this r e spect, it is interesting to note that 
INoordermeer & VerheiienI (l2007l) found that, depending on 
which velocity measure they used, there is a subset of massive 
spiral galaxies rotating ’too fast’ (or equivalently, that are ’too 
dim’), causing a ’kink’ or a change of slope in the relation. This 
effect is smallest for their sample if they use the circular velocity 
derived for large radius as opposed to one based on the maximum 


Fig- 9. Res idual magnitudes with respect to the CO TFR from 
iDavis et H] ( 1201 Ih as function of M/Ljam- The dashed line indicates 
how a 20% distance error moves an object in the plot. Such an error 
corresponds to an error in absolute magnitude of about 0.4 mag and of 
0.07 in log MIL,AM 


velocity more representative for the inner regions. They interpret 
this as an indication that the shape of the rotation curve changes 
systematically with mass in the sense that more massive galaxies 
tend to have a rotation curve that is slightly declining while it is 
more flat for less massive spirals. Above we found evidence that 
the rotation curves of our early-type galaxies are behaving simi¬ 
larly to the massive spiral galaxies of iNoordermeer & VerheiienI 

(l2(m . __ 

Figure|2]shows the lNoordermeer & VerheiienI d2007h data for 
spiral galaxies, based on the asymptotic velocity, compared to 
our data for the ATLAS^*^ galaxies. Our galaxies still show an 
offset to fainter magnitudes despite the fact that we also use the 
circular velocity at large radius so the effects of decreasing cir¬ 
cular velocities should be largely taken out. This suggests that 
differences in stellar population are the main cause for the ob¬ 
served offsets. 

To further investigate this, we note that Fig. [T] as well as 
the results given in Table 2, show that the slope of the TFR for 
ETGs is flatter than that of later-type galaxies, pointing to some 
systematic effect. FigureISjshows the offset o f our galaxies from 
the TFR of iNoordermeer & Verheiieij[j 2007h as a function of the 
r-band mass-to-light ratio derived bv ICappellari et akl (1201 3l) us¬ 
ing dynamical modelling of the observed stellar kinematics (see 
below). There is a clear trend visible in this figure, which un¬ 
derlines that differences in stellar populations play a major role 
when comparing the TFRs of early-type and late-type galaxies 
and that a more accurate scaling relation, based on mass instead 
of light, exists. 

This is also suggested by Fig. 0 where we plot the residual 
(in magnitude) of our galaxies from the CO TFR as function of 
the MjL as in Fig.|3 It is clear that also here a systematic trend 
is present, albeit with a d ifferent slope than in Fig. | 8 ] due to the 
fact that the CO TFR of IDavis et al. ( 20111) has a flatter slope 
[Nc 


than the TFR given by INoordermeer & V^heiienl (l2007h . The 
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Fig. 10. The baryonic TFR derived using dynamically determined M/L’s (left) and Ml Vs derived from modelling the star formation history (right). 


fact that a correlation exists between the offset from the TFR 
and the MjL of the stellar population suggests that the scaling 
relation between baryonic mass and circular velocity will have 
smaller scatter than the one involving light and circular velocity. 

Before we investigate this, one should note, however, that 
both the absolute magnitude as well as M/Ljam are distant de¬ 
pendent quantities and therefore both suffer from distance er¬ 
rors which could, in principle, cause a correlation of the kind 
seen in Fig. 0 In fact, the absolute magnitude depends on 
while M/Ljam depends on D In Fig. |9]is indicated, with the 
dashed line, how distance errors would move data points through 
the hgure. Figure 0 shows that the data points follow a similar 
slope as would be expected from distance errors, but that the ob¬ 
served spread is much larger than expected from the errors in 
distance for our sample galaxies. If the observed spread would 
be entirely due to distance errors, it would imply that the dis¬ 
tances are uncertain up to a factor 2 which is much larger than 
the typical ~20% dis tance error estimated in Sec. |2] (see also 
ICapDellari et al.ll201 Ih . 


5. The baryonic TFR of early-type galaxies 


The TFR is an empirical relation between circular velocity 
and luminosity which may result from a more fundamental 
relation between galaxy properties. In this context, a num¬ 
ber of authors have investigated the existence of a bary¬ 
onic TFR^_where_total_baryonicjim^ replaces luminosity (see 
e.g. iMcGaugh et al.l l2000t iMcGaughI l2()12l) . It has been sug¬ 
gested that the baryonic TFR is a more fundamental rela¬ 
tion since it appe ars to have less sca tt er than the classi¬ 
cal TFR (see e.g. iGuroyich et al. 20041: Begum etakl l2008l 


Noordermeer& VerheiienI l2007l IZaritskv et^ 20141) . Almost 
all earlier studies of the baryonic TFR have focussed on spiral 
galaxies. If, as we concluded in the previous Section, systematic 
differences in stellar population are important for explaining sys¬ 
tematic trends in the residuals of ETGs from the standard TFR, 
it is interesting to include early-type galaxies when considering 
the baryonic TFR. 


We compute the baryonic mass of a galaxy as the sum of 
the cold-gas mass and the stellar mass, although for our sample 
galaxies, the stellar mass is the dominant contribution; e.g., the 
H i-to-stellar mass ratio ranges between 0.1 to 8 % (2.2 % on 
av erage). The atomic gas mass is computed as the H i mass given 
in ISerra et ^ (l2012l) . multiplied by a factor of 1.4 to take into 
account the presence of helium and metal s. Moreover we ad d 
the mass of the molecular component from lYoung et akl (1201 ll) . 
which is also only a very small contribution. 

S tellar masses are calcu lated using the r-band luminosities 
from ICappellari et"^ (1201 3l) and multiplying them with a mass- 
to-light ratio MjL;. We consider two different estimates of this 
ratio: hrstly, the mass-to-light ratios resulting from JAM mod¬ 
elling of t he stellar kinema t ics, de rived for the full ATLAS 
sample by ICappellari et al.l (1201 3l) where we use their results 
from the self-consistent JAM models (their model A). The sec¬ 
ond set of mass-to-light ratios we use are based on the star- 
formation histories (SFHs) derived using a Salpeter IMF (Cap- 
pellari et al. 2013a). Note that Cappellari et al. (2012), based on a 
comparison of these sets of mass-to-light ratios, found that there 
are systematic variations in the IMF among early-type galaxies 
and that assuming a Saltpeter IMF for all galaxies will lead to 
systematic effects in the estimates of the galaxy masses. 

Figure [To] shows the baryonic TFRs derived for the two dif¬ 
ferent sets of mass-to-light ratios, and Table 2 shows the param¬ 
eters of the best-fitting relations for 3 different sets of assump¬ 
tions.. The baryonic relation derived using the dynamically de¬ 
termined M/Ljam shows less scatter than the one based on the 
MjL derived from modelling the star formation histories of the 
galaxies. In addition, the scatter in the SFH based baryonic TFR 
is very similar (when converted to optical magnitudes) to that 
of the standard TFR while the scatter in the JAM-based TFR is 
signihcantly smaller. Both these features are exactly what is ex¬ 
pected if the underlying baryonic TFR is tight while the main 
sources of scatter are the uncertainties in the distances to the 
galaxies. Given that M/Lsfh is distance independent, whereas 
M/Ljam has a D ' distance dependence, the errors on the masses 
using M/Lsfh vary as whereas those on the masses using 
M/Ljam vary only linearly with D. Therefore a 20% distance 
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Fig. 11. Baryonic TFR where the galaxy m asses for our sample galaxi e s are estimated using M/Ljam- The red symbols re present our ATLA S^° 
sample, the black symbols are the data fr omlNoo rdermeer & Verheiienl (l2007h. while the open symbols show the data from lMcGaughl (1201 2h for 
gas-dominated galaxies. For the d ata from lNoordermeer & Verheiien l l2007h in the left-hand panel MILk - 0.8 MqILq was used, as in the original 
iNoordermeer & Verheiienl ( l2007ll paper, while in the right-hand panel M/L^- = 0.54 Mq/Lq was used. The scaling of the lefthand figure is the 
same as that of Fig.[7]to facilitate easy comparison. 


error implies a typical uncertainty of ~0.17 dex on logMbar.SFH 
(or of 0.43 when converted to magnitudes, as in Table |2]l and 
of ~0.09 dex on log Mbar.jAM (or 0.23 when converted to magni¬ 
tudes). The observed scatter of the baryonic TFRs of our sample 
is only slightly larger than these estimates. Given that the galaxy 
masses based on M/Ljam suffer less from distance errors, and 
also that they suffer less from systematic errors than masses de¬ 
rived using M/Lsfh, in the following we only further consider 
the baryonic masses using M/Ljam- 

Table|2]gives the results of fits, of the form 

logMbar = a (log W - 2.6) -F b. 


0.43 mag (compared to 0.62 mag for the standard TFR). Re¬ 
cent work has su ggested that the disks of spiral galaxies are 
not maximal (e.g., Bershadv et al.ll2010l : lMartinsson et alJl2013l : 


iMcGaugh & Schomberl)l2oT5 ) so this offset may simply reflect 

that using MILk - 0.8 MqILq overestimates the stellar masses 
of spiral galaxies. T ransforming the offset between ou r baryonic 
TFR and the one of iNoordermeer & Verheijet^ (l2007h into a re¬ 
duction of MILk suggests that one should lower this quantity 
with a factor 0.67 so that MILk - 0.54 MqILq would be a more 
appropriate value for sp iral galaxies, in line with the results of 
some recent work (e.g. Bershadv et al.l 1^1 Ot iMartinsson et al.l 
120 131: IMcGaugh & Schombertll20L5h. 


to our data for a number of cases. This table shows that, if both 
the slope a and the intercept b both are free parameters, the slope 
is flatter than that of the baryonic TFR for spiral galaxies der ived 
by Noordermeer & Verheiienl (l2007h and iMcGaughI (1201 2h . al¬ 
though, given the small sample size and the limited range in log 
W covered by our sample, this result is somewhat uncertain. 

In Fi g, fin we compare the baryonic TFR with that de¬ 
rived by (Noordermeer & Verheiienl (l2007h for their sample of 
massive spiral galaxies. In deriving baryonic masses, they as¬ 
sumed, in a quite different approach to what we have chosen, 
a fixed MILk - 0.8 MqILq for all galaxies in their sam¬ 
ple, a value motivated by the maximum-disk scenario. Despite 
the different approaches, a comparison of Fig. [TT] with Fig. [T] 


shows that our earl y-type galaxies lie much closer t o the spi¬ 
ral baryonic TFR of iNoorderm eer & Verheii e^(l2007h than the y 
do to the standard TFR of INoordermeer & Verheiienl (l2007ll . 
W hen we fit our baryonic TFRs wi th the same slope as found 
by Noordermeer & Verheiienl (l2007l) for their relation, we And 
an intercept of 10.71 + 0 .03 co mpared to 10.88 + 0.02 for 
Noordermeer & Verheiienl (l2007ll sample. This implies an off¬ 
set between the baryonic TFR of spirals and of ETGs of about 
0.17 + 0.04 dex in mass, or when expressed as magnitudes, it is 


In the right-hand panel of Fig. [TT] we explore how our re¬ 
sults compare to the baryonic TFR for low-mass, gas-rich galax- 
ies. We compare our data with the baryonic TFR presented by 
lMcGaughl(l2012l) . This latter sample consists of particularly gas- 
rich galaxies (Mgas > Mjtar) to minimise the influence of errors 
in the assum ed the stellar ma ss-to-light ratios which haye been 
estimated by lMcGaughl (1201 2h using a Kroupa initial mass func¬ 
tion. ]pJhej;ighFhand_panel_ofFig;^T] we also include the data 
from lNoordermeer & Verheiienl (l200% where we now assumed a 
Axed MILk - 0.54 MqILq for all their galaxies. FigurefTTIshows 
that, despite the different methods used for computing stellar 
masses, the combined dataset produces a well-defined baryonic 
TFR that seems to be yalid for small gas-rich galaxies, massiye 
spiral galaxies as well as early-type galaxies. iMcGaughI (1201 2h 
fitted a relation of the form log Mbar = 3.82 (log W-2.6 )-f 10.79. 
When we fit a baryonic TFR with the same slope to our data 
(Table |2ll we find a small, but not yery significant offset, of 
0.06 + 0.06 dex for our galaxies. This is quite different from 
the situation for the standard TFR where the early-type galaxies 
lie clearly offset from the spiral TFR. The well-defined baryonic 
TFR shows that the offsets from the standard TFR are due to 
different stellar populations between the two classes of galaxies 
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and that the fundamental relation is more likely to be the bary- 
onic TFR. 


6. Conclusions 

We study the H i /T-band TFR and the baryonic TFR of a sam¬ 
ple of 16 ETGs taken from the ATLAS^*^ sample which all have 
very regular H i disks that extend well beyond the optical body 
5Refr)- We use the kinematics of these disks to estimate the 
circular velocity at large radius for these galaxies. For galaxies 
that are sufficiently well resolved spatially, we use the veloc¬ 
ity fields, including a harmonic analysis, to find the best-fitting 
t^circ for the outermost point of the rotation curve. In three other 
cases, we model the full data cube in order to derive the kine- 
matical information of the galaxy, and in two additional cases, 
we use kinematical data from the literature. We use these circu¬ 
lar velocities to construct the traditional fT-band TFR as well as 
the baryonic TFR. In the case of the baryonic TFR, we use both 
mass-to-light ratios derived from stellar kinematics and from star 
formation histories. We find the following; 


1. Comparison with published TFRs derived for samples of spi¬ 
ral galaxies suggests that the TFR for ETGs is offset by about 
0.5-0.7 mag in the sense that early-type galaxies are dimmer 
for a given circular velocity. The residuals from the standard 
TER correlate with estimates of the M/L of the stellar popu¬ 
lations, suggesting that the offset is mainly driven by differ¬ 
ences in stellar populations. 

2. There is a small offset between the TER derived from circular 
velocities at large radius from our H i data with the relation 
derived for the ATLAS^° sample using CO data covering the 
galaxies’ inner regions (< IReft)- This offset suggests that 
the circular velocities at large radius are about 10% lower 
than the CO circular velocities near 1 Reft- Such a decrease 
is consistent with rece nt determinations of th e shape of the 
mass profiles of ETGs (ICappellari et al.ll2015l) . 

3. The baryonic TER is distinctly tighter than the standard TER, 
in particular when using the mass-to-light ratios based on 
dynamical models. In both cases the scatter is mainly due 
to distance uncertainties. Erom t he offset we find between 
our ba ryonic TER with that of iNoordermeer & Verheijeni 
(l2007h derived for spiral galaxies, we find that the ETGs fall 
on the spiral baryonic TER if one assumes M/Lk - 0.54 
MqILq for the stellar populations of the spirals. Such a value 
would imply that the disks of spiral galaxies are about 60- 
70% of maximal, as is also found by other recent studies of 
the dynamics of spiral galaxies fe.g. iBershadv et al. 


iMartinsson et aljr2013l; iMcGaugh & Schombertll20'l5[) . The 
galaxies of our sample also fall on the baryonic TER de¬ 
rived for gas-dominated, late-type, lower-mass galaxies. The 
fact that there is no significant offset between the baryonic 
TERs of spirals and ETGs suggests that the differences in 
concentration of the mass distributions of these two types of 
galaxies have little effect on the baryonic TER. Our analy¬ 
sis increases the range of galaxy morphologies for which the 
baryonic TER holds, strengthening previous claims that this 
is a more fundamental scaling relation than the classical re¬ 
lation. 
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